Universal property of the information entropy in atoms, nuclei and atomic clusters 



S.E. Massen and CP. Panos 
Department of Theoretical Physics Aristotle University of Thessaloniki, GR 54006 Thessaloniki, Greece 



oo 



> 

o 

(N 
O 
OO 
ON 



X 



The position- and momentum-space information entropies 
of the electron distributions of atomic clusters are calculated 
using a Woods-Saxon single particle potential. The same en- 
tropies are also calculated for nuclear distributions according 
to the Skyrme parametrization of the nuclear mean field. It 
turns out that a similar functional form 5 = a + fe In A'^ for the 
entropy as function of the number of particles A'^ holds ap- 
proximately for atoms, nuclei and atomic clusters. It is con- 
jectured that this is a universal property of a many-fermion 
system in a mean field. It is also seen that there is an anal- 
ogy of our expression for S to Boltzmann's thermodynamic 
entropy S — klnW . 

PACs: 89.70.+C; 36.40.+d; 31.10.+z; 21.60.-n 

Information-theoretical methods have played in recent 
years an important role in the study of quantum mechan- 
ical systems in two cases: first in the clarification 
of fundamental concepts of quantum mechanics and sec- 
ond in the synthesis of probability densities in position 
and momentum spaces. In the first case an important 
step was the discovery of an entropic uncertainty rela- 
tion (EUR) by Bialynicki-Birula and Mycielski ^ which 
for a three-dimensional system has the form: 



(1) 



(see also ref . for the one-dimensional case) . In (|l|) Sr 
is the Shannon information entropy in position-space: 



Sr = - J pir)\np{r) df (2) 
Sk is the corresponding entropy in momentum-space: 

Sk = - [ n{k)\nn{k) dk (3) 



and p{r)^ n{k) are the position- and momentum-space 
density distributions respectively, which are normalized 
to one. However, for a normalization to the number of 
particles N, the following EUR holds Q: 

Sr + Sk> iN{l + In tt) - 2iV In = 6.434iV - 2A^ In N 

(4) 

Inequality ( 1 ) , for the information entropy sum in con- 
jugate spaces, is a joint measure of uncertainty of a 
quantum mechanical distribution, since a highly local- 
ized p{r) is associated with a diffuse n(fc), leadiiig to 
low Sr and high 5*^ and vice-versa. Expression (|l|) is 
an information-theoretical uncertainty relation stronger 
than Heisenberg's We also note that expression (||) 



does not depend on the unit of length in measuring p{r) 
and n{k) i.e. the sum Sr + 5*^ is invariant to uniform 
scaling of coordinates. 

Gadrc |^ derived the following approximate expression 
for the information entropies of electron distributions in 
atoms: 



Sr + Sk ~QmN~N\nN 



(5) 



using Thomas-Fermi theory and Gadre et al |g| derived 



Sr + Sk^ Q.2b7N - 0.993A^lniV 



(6) 



with Hartee-Fock calculations. Here, N is the number of 
electrons. 

Panos and Massen Q found the following expression 
for nuclear distributions, employing the simple harmonic 
oscillator (HO) model of the nucleus: 



Sr + Sk^ 5.287iV - l.lSiVlnTV 



(7) 



where N is the number of nucleons in nuclei. Relations of 
the same functional form hold for Sr and Sk separately 
but the important quantity is Sr -\- Sk- 

There is a striking similarity of (5), (6) and (7) with 
the EUR (4), indicating that the functional form 



S ^ aN + hN\nN 



(8) 



is universal for a many-fermion system in a mean field. 

However, the above relations were derived for a nor- 
malization of p{r) and n{k) to the number of particles 
N . In the following we find it more convenient to nor- 
malize to one . There is a simple relationship between 
the two cases and we can easily transform one case to the 
other according to the relations: 

r ,1 Sr\norm = N] , 

Srlnorm = 1 = — +\nN 

^ ' N 



^ r 1 Sklnorm — N] , ,^ 
Sk [norm = 1] = — — + \nN 

Hence, we have for normalization to one, the following 
expressions: 

Sr + Sk — 6.65 -f In A'^ (atoms, Thomas — Fermi) (9) 



Sr + Sk^ 6.257 + 1.007 IniV (atoms, Hartee - Fock) 

(10) 

5r + S'fc - 5.287 + 0.870 In iV (nuclei - H.O.) (11) 
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In the present letter we extend our calculations for 
two other cases: the distribution of the valence electrons 
in atomic clusters using a Woods-Saxon single particle 
potential and the nuclear distribution in nuclei employing 
the Skyrme parametrization of the nuclear mean field. 

In atomic (metallic) clusters the effective radial elec- 
tronic potential was derived by Ekardt in his 
spherical-jellium-background-model study of the self- 
consistent charge density and the self-consistent effective 
one-particle potential, using the local density approxi- 
mations. Ekardt's potentials for neutral sodium clusters 
were parametrized in ref. jl^ by a Woods-Saxon poten- 
tial of the form: 



Vws{r) 



Vo 



1 + exp[--^] 



(12) 



with Vo ^ 6 eV, R = tqN^^^, ro = 2.25 1 and a = 
0.74 A. For a detailed study regarding the parametriza- 
tion of Ekardt's potentials see ref. |]l3|| . 

We solved numerically the Schrodinger equation for 
atomic clusters with Z = 8, 18, 20, 34, 40, 58, 68, 70 
valence electrons in the potential ( p^ and found the 
wave functions of the single-particle states in configu- 
ration space and by Fourier transform the corresponding 
ones in momentum space. Using the above wave func- 
tions, we calculated the electron density p{r) in position 
space and n{k) in momentum space, which were inserted 
into equations (||) , (||) and gave us the values of the in- 
formation entropies Sr and Sk ■ Then we fitted the form 
5' = a + 6 In to these values and obtained the expres- 
sions: 



Sr ^ 4.133 -f 0.934 IniV 
Sk 1.563 - 0.027 In iV 
S*^ + S'fc~ 5.695 + 0.907 In iV 



(13) 
(14) 
(15) 



Next the nuclear densities p(r) and n{k) for several nu- 
clei were obtained with Hartee-Fock calculations using 
the Skyrme parametrization of the nuclear mean field. 
There are various parametrizations of the Skyrme inter- 
action, but they affect slightly the information entropies 
|||. Thus we used the SKIII interaction Q. Finally, we 
fitted the form 5 = a + 5 In iV to the values obtained with 
SKIII interaction and we found the expressions: 



Sf 



Sr ~ 3.395-^0.767 In TV 
Sk 1.929 -I- 0.091 In iV 
Sk 5.319 + 0.860 In iV 



(16) 
(17) 
(18) 



The fit is in reasonably good agreement with its H.O. 
counterpart (comparison of relation (18) to ([Tl])), though 
the individual entropies Sr and Sk do not match with the 
respective HO ones § that well. It seems that there is a 
delicate balance between the coordinate and momentum 
spaces, so that the interesting quantity is the sum S = 
Sr + Sk (the net information content of the system) and 
not the individual entropies Sr and Sk- 




FIG. 1. Information entropy S as function of the number of 
particles according to our fitted expression S = a + bin N 
for atoms, clusters and nuclei. The corresponding values of 
our numerical calculations are denoted by solid circles (clus- 
ters) and open circles (nuclei). 

In figure 1 we plot our fitted form S = Sr + Sk = 
a + 6 In TV for atoms (with Hartee-Fock, relation (IC), 
upper curve), atomic neutral Na clusters (relation (15), 
middle curve) and nuclei (with Skyrme, relation [li), 
lower curve). These lines correspond to our fitted ex- 
pressions, while the corresponding values of our numeri- 
cal calculations are denoted by solid circles (clusters) and 
open circles (nuclei with SKIII interaction). 

Concluding, in the present letter we derive an in- 
teresting characteristic of information entropies Sr and 
Sk for various systems i.e. atoms (Thomas-Fermi the- 
ory, Hartee-Fock), nuclei (Harmonic Oscillator model, 
Skyrme force) and atomic clusters (Woods-Saxon poten- 
tial) . For all of these systems the entropies can be repre- 
sented well by a function, which incorporates IniV lin- 
early i.e. S = a + bliiN where N is the number of 
electrons in atoms or nucleons in nuclei or electrons in 
atomic clusters. We may conjecture that this is a univer- 
sal property of a many-fcrmion system in a mean field. 
As stated in ref. the information entropies seem to 
be a hidden treasure, as yet remains mostly unexplored. 

A final comment seems appropriate: There is an anal- 
ogy of our expression S = a + bin N to Boltzmann's 
thermodynamic entropy 5 = fc In (fc is Boltzmann's 
constant). This can be seen as follows: Our relation can 
be written in the form: 

5 = 61n(iV/A^o) (19) 

where I/Nq — e""!^ . The above entropy is the information 
entropy S"™/, which is related to the physical entropy 



Sphya through Jaynes' relation: S^hys 
case, we have: 

S^hy, = fc51n(iV/A^o) 



kS, 



inf ■ 



In our 



(20) 



Here b and Nq depend on the system under consideration 
(atom, cluster or nucleus), but this dependence is not 
strong. One can take that 6 « 1 (actually equals 1.007 for 
atoms, 0.907 for clusters and 0.860 for nuclei) and I/Nq « 
500 (because 1/iVo = 500 for atoms, 485 for nuclei and 
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533 for clusters). Thus we obtain the following rough 
estimate for the entropy of a many-fermion system: 

Sphys ^ k\n{500N) (21) 

The number 500A^ corresponds to the number of 
equiprobable microstates W of Boltzmann's relation. 

Hence, starting from Shannon's information entropies 
of quantum-mechanical distributions, we arrived at re- 
lation which can be considered as a quantum- 
mechanical analogue to Boltzmann's entropy 5 = fc In W. 
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